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ABSTRACT 

It is generally accepted that the pressure-dilatation, which is an additional compressibility 
term in the turbulence transport equations, may be important for high-speed flows. Recent 
direct simulations of homogeneous shear turbulence have given concrete evidence that the 
pressure-dilatation is important insofar that it contributes to the reduced growth of turbulent 
kinetic energy due to compressibility effects. The present work addresses the problem of 
modeling the pressure-dilatation. We first isolate a component of the pressure-dilatation 
which exhibits temporal oscillations and, using direct numerical simulations of homogeneous 
shear turbulence and isotropic turbulence, show that it has a negligible contribution to 
the evolution of turbulent kinetic energy. Then, an analysis for the case of homogeneous 
turbulence is performed to obtain a model for the non-oscillatory pressure-dilatation. This 
model algebraically relates the pressure-dilatation to quantities traditionally obtained in 
incompressible turbulence closures. The model is validated by direct comparison with the 
pressure-dilatation data obtained from the simulations. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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1 Introduction 


The pressure-dilatation appears as an explicit compressibility term in the equations for mean 
temperature and turbulent kinetic energy. It has been generally recognized that this term 
may be large in high-speed flows and therefore requires consideration in compressible turbu- 
lence closures. In Sarkar, Erlebacher, Hussaini and Kreiss 1 , it was found that the pressure- 
dilatation was smaller than the compressible dissipation in direct simulations of isotropic 
compressible turbulence which suggested that the presure-dilatation can be absorbed into 
the model for the compressible dissipation derived therein. However, subsequent direct sim- 
ulations of homogeneous shear flow by Blaisdell et al . 2 and Sarkar, Erlebacher and Hussaini 3 
showed that the pressure-dilatation is comparable to the compressible dissipation and con- 
tributes to the reduced growth of turbulent kinetic energy induced by compressibility. This 
has motivated a revisit to the issue of modeling the pressure- dilatation in the present paper. 
Recently, Taulbee and Van Os cl ol 1 and Zeman 5 have also considered this problem. Taulbee 
and VanOsdol related the sum of the pressure-dilatation and compressible dissipation to a 
model involving the density variance and the divergence of the mean velocity. A separate 
modeled transport equation was proposed for the density variance. Zeman 5 equated the 
pressure-dilatation to the time rate of change of pressure-variance and then, extending the 
work by Sarkar et al . 3 on equilibration of the compressible pressure variance on the fast 
acoustic time scale, proposed a transport equation for the pressure variance. The objective 
of the present work is to deduce a model for the pressure-dilatation which is algebraically 
related to quantities obtained by incompressible turbulence closures, in contrast to the pre- 
vious models, which require new transport equations. 

We consider the behavior of the pressure-dilatation in homogeneous flows. Direct simu- 
lation results are presented and compared for two flows - homogeneous shear turbulence and 
decaying isotropic turbulence. The DNS results show that the contribution of the pressure- 
dilatation to the turbulent kinetic energy evolution is more important in the case of shear 
turbulence than in decaying isotropic turbulence. We also find the somewhat surprising re- 
sult that the major contributor to the pressure-dilatation is not the compressible pressure, 
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but is the incompressible pressure associated with the solenoidal velocity. The governing 
equations are analyzed for the case of homogeneous flow to obtain a formal expression for 
the pressure-dilatation from which a model is deduced using scaling arguments. 

2 DNS results 

A spectral collocation method along with a third-order accurate Runge-Kutta time advance- 
ment was used to solve the compressible Navier-Stokes equations. Details of the calculations 
are available in Sarkar, Erlebacher and Hussaini 3 for homogeneous shear flow, and in Er- 
lebacher, Hussaini, Kreiss and Sarkar 6 for decaying isotropic turbulence. A uniform 96 3 
grid was used for discretizing the flow domain. Table 1 shows some of the parameters of the 
three shear flow cases (Sl,S2 and S3), and three decaying isotropic cases (D1,D2, and D3). In 
Tables 1 and 2, S and v denote the shear rate and kinematic viscosity used in the (nondimen- 
sional) compressible Navier-Stokes equations, while Rx,o a-nd Mt , o denote the initial values 
of the Taylor microscale Reynolds number and turbulent Mach number respectively. Note 
that Re x = qX/v where q = and A = q/\Ju/M , while M t = qfc where c = \J-yRT is 

the mean speed of sound. In order to minimize the introduction of compressible effects due 
to initial conditions, all the cases start with incompressible data, that is, the velocity field 
is divergence-free ( d 1 — V.u / = 0), the pressure field satisfies the usual Poisson equation for 
incompressible flows, and density fluctuations p' = 0. The initial temperature is calculated 
from the equation of state using the known pressure and density fields. 

The homogeneous flows considered here have temporally evolving turbulence statistics. 
Fig. 1 shows the evolution of the pressure-dilatation p'd! as a function of normalized time 
St for the two homogeneous shear cases SI and S2. Two trends in the behavior of p'd 1 
are evident from Figs, la-b; first, the pressure-dilatation develops pronounced oscillations 
in time and second, it is more negative than positive. These trends are common to all the 
simulations of homogeneous shear flow that we have performed. Such oscillations in statisti- 
cal quantities are neither expected nor encountered in incompressible flows and complicate 
issues in compressible turbulence modeling. By numerical experiments, it was found that 
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the nominal time period of the oscillations decreased approximately linearly with the speed 
of sound. This suggested that one could isolate the oscillatory part of p'd' by decomposing 
the fluctuating pressure p into the sum of an incompressible part p ^ and a compressible 
part p c . The component p 1 ' is associated with the incompressible velocity field u 7 which is 
divergence-free (V.u 7 = 0) and satisfies the usual Poisson equation 


vy ' = 


^ P U m,n U n,m P U m,n U n,rt 


( 1 ) 


and the remainder p' — p 1 ' is the compressible pressure p c ' . Since 


P = P 1 ' + P C ' 


we have 


j/d* = p I! d f + d f 


Fig. 2a shows the evolution of p c 'd' and p 7 'd' for case Si. The oscillations are substantial 
only for p c ‘d', and furthermore, the peaks and valleys in the evolution of p c ‘ d' in Fig. 2a 
seem to be much more symmetric around the origin than those in p'd 7 in Fig. la. 

One of the important effects of p'd' is its influence on the budget of the turbulent kinetic 
energy K . The equation for K in homogeneous turbulence is 


P dt^ K ^ = ~pP - pt a - pe c + p'd' 


( 2 ) 


where V = -Uiju'/u'j is the production, e, the solenoidal dissipation, e c the compressible 
dissipation, and p'd' the pressure-dilatation. The overbar denotes a conventional Reynolds 
average, while the overtilde denotes a Favre average. A single superscript ' represents fluctua- 
tions with respect to the Reynolds average, while a double superscript " signifies fluctuations 
with respect to the Favre average. In order to gauge the relative importance of the two com- 
ponents p c 'd' and p 7 'd' of the pressure-dilatation in the evolution of the turbulent kinetic 
energy, we plot the integrals J^d'dt and / p c 'd'dt in Fig. 2b. Fig. 2b shows that the inte- 
grated contribution of p 7 d' is about an order of magnitude larger than that of p^' d' . Thus, 
even though the extrema of p c 'd' are larger than p 7 'd', the time-integrated contribution of 
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these oscillations in p c 'd’ is negligible compared to the contribution of p r d'. Fig. 3 and 
Fig. 4, which contrast the behavior of ^ and for two other homogeneous shear flow 

cases S2 and S3, confirm the aforementioned trends observed in case SI. (See Table 1 for 
the initial conditions of cases S2 and S3). It seems that for homogeneous shear flow, p v d! 
dominates p^d/ in the kinetic energy equation. It is interesting that even though both the 
dilatation and the compressible pressure are associated with the same hyperbolic system of 
equations governing the compressible mode, the correlation of the dilatation with the com- 
pressible pressure plays a smaller overall role than the correlation of the dilatation with the 


incompressible pressure. 

Figs. 5 and 6 show DNS results on the pressure-dilatation from the decaying isotropic 
turbulence cases D1 and D2 (see Table 2 for initial conditions). Both p v d! and p c ' d! in 
Fig. 5a and 6a show sharp transients initially whereby the pressure field is appropriately 
redistributed into incompressible and compressible parts. The integrated p 1 ' d' is much larger 

than the integrated W* in Figs. 5b and 6b J ust as in the bom °g eneous shear flow case ' In 
Figs. 5a and 6a, the term p c 'd' becomes approximately zero after a short initial transient. 
This is in agreement with the exact solution derived in Sarkar et al. 1 for the linearized 
equations for the compressible mode which predicted that p c> d' -> 0 after a transient on 
the acoustic time scale. The DNS reveal an important difference in the behavior of j id! 


in isotropic turbulence with respect to homogeneous shear. The term p'd' is predominantly 
positive in the case of decaying isotropic turbulence, in contrast to its predominantly negative 
behavior in homogeneous shear turbulence. The different signs of the pressure-dilatation have 
been explained by a theoretical consideration of the equations of the pressure variance and 


density variance by Sarkar, Erlebacher and Hussaini 3 . 

To summarize, the DNS results show that the pressure-dilatation is predominantly neg- 
ative and has pronounced oscillations in homogeneous shear flow, while it is predominantly 
positive in decaying isotropic turbulence. We find that the oscillations in Jd! are confined 
to the component ~f r d' of the pressure-dilatation associated with the compressible pressure, 
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and furthermore that, because of self-cancellation, the contribution of this oscillatory com- 
ponent to the development of the kinetic energy is negligible. Therefore, it seems that only 
the component p 1 ' d' of the pressure-dilatation requires modeling, which we proceed to do in 
the next section. 


3 Modeling the pressure-dilatation 


Consider the Poisson equation Eq. (1) for the incompressible pressure. As is often done for 
pressure-strain modeling (see Lumley 7 , Reynolds 8 ) for incompressible flows, the incompress- 
ible pressure can be split into a rapid part p & which reacts instantaneously to a change in 
the mean velocity gradient and a slow part p s \ Thus 


P '' = ? k 


+p 


s' 


( 3 ) 


where 


and 


V 2 p*' = -2 


v 2 / = -w m y n , m ( 5 ) 

Eq. (4) can be exactly solved by Fourier transforms for homogeneous turbulence to give 


( 4 ) 


P R ' = 


( 6 ) 


Note that <j) denotes the Fourier transform of . From Eq. (6) it follows that 

puL,J ^E^fdk (7) 


Here E^f denotes the spectrum of the mixed Reynolds stress tensor u R uf. Eq. (7) is an 
exact expression for the rapid pressure-dilatation p R ' d.' and, though cumbersome, can be 
used for obtaining a simple model. By inspection of Eq. (7), it is clear that because of 
the dependence on the local energy spectrum tensors, a transport equation is required for 
a general representation of However, we will attempt to obtain an algebraic rather 
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than a differential model for compressibility terms. This is possible if there are equilibrium 
scalings in the flow, and the ensuing model will be useful if the additional compressibility 
correlations do not dominate the incompressible terms in a given transport equation. 

Define the tensor A mn by 

<4m» = / ^Egdk (8) 

which enables the following compact representation of Eq. (7) of the rapid pressure-dilatation 

W* = 2 ( 9 ) 

The simplest dimensionally consistent form for A mn which has the correct dependence 
on uj' and uf’ is 

( 10 ) 

where <x 2 is a dimensionless parameter which in general is a function of the actual shape 
of the energy spectrum tensors for the incompressible and compressible velocity. We have 
found through analysis and DNS that for homogeneous shear turbulence 

u c = 0(M t )u I (11) 

where u G and u 1 respectively denote the L 2 norms of the compressible and incompressible 
velocity fields. In Sarkar et al.^, we had shown that the dilatational velocity has a fast time 
scale which is 0(M t ) times the solenoidal velocity’s time scale, and therefore the correlation 
between u l m ' and u G> should be prorated by a factor of M t . Using Eq. (11) and prorating 
the mixed correlation, we obtain from Eq. (10), 

Amn = a2Mt U m U n ( 12 ) 

Substituting Eq. (12) into Eq. (9) gives 

p^d’ = -a 2 M?V (13) 

where V is the production Eq- ( i3 ) ma y be used as a model for the rapid ' 

pressure dilatation when using a conventional Reynolds- averaged system of equations. We 
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prefer using a combination of Favre and Reynolds averages (eg. see Sarkar and Balakrishnan 9 
for compressible flow calculations, and for such a system of equations we propose using 
Eq. (13) with the production being V = Since p' j~p — 0(Af 2 ), changing from the 

Reynolds-averaged definition to the Favre-averaged definition of the production V introduces 
terms of higher order in which may be neglected. 

Let us now consider the remaining part of the pressure-dilatation, the slow pressure- 
dilatation p s 'd'. After using Fourier transforms to solve Eq. (5) for the slow pressure, we 
obtain the following expression for the slow pressure-dilatation, 

p s ’d' --p j (iuiu/u^* - iuplj 1 *^) dk (14) 

Here <^* denotes the complex conjugate of the Fourier transform <^. An order of magnitude 
analysis of the r.h.s. of Eq. (14) gives 

?V = T>0(^)0(^-)M, (15) 

In Eq. (15), / denotes the integral length scale of the turbulence, and the last factor M t arises 
from the disparity between the scales of the incompressible and compressible fields. Using 
Eq. (11) for the scaling of u c , and noting that e s = 0(u /3 //), we obtain the model 

p s 'd f = a 3 pc 3 M? (16) 

Combining Eqs.(13) and (16), we have the following model for the incompressible pressure- 
dilatation 

p v d f = -a 2 pVMt + a 3 pe s M? (17) 

We will now use the DNS data to verify the functional dependence stated in Eq. (17) and 
also to calibrate the model coefficients a 2 and a 3 . 

Because the production V = 0 in decaying isotropic turbulence, the variation of the 
incompressible pressure-dilatation with e 9 can be verified using DNS of isotropic turbulence. 
The ratio p p d f /(’pe 3 M^) is shown as a function of non-dimensional time in Fig. 7a. This 
ratio reaches an equilibrium value by a time of 0.25, substantiating the validity of the second 
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term in Eq. (17). Based on the DNS value of the equilibrium ratio, the model coefficient q 3 
in Eq. (17) is taken to be 0.2. The remaining part of the model for the pressure-dilatation 
is calibrated against simulations of homogeneous shear flow. Fig. 7b shows that, in accord 
with our model, the rapid part of the pressure-dilatation scales as ~pPMf. The ratio (p 1 ' d! — 
0.2 e,M t 2 )/(pT > M i 2 ) reaches an approximate equilibrium value of —0.4, suggesting that the 
model coefficient a 2 = 0.4. 

Finally, after using the arguments of section 2 to neglect p c> d! relative to p 1 ' d' , the model 
for the pressure-dilatation becomes 

p'd' = -a 2 pVMt + a 3 pe a M j (18) 

where V = —uiju'-u" is the production of kinetic energy, M t = \J2kJ^bT the turbulent 
Mach number, e s the solenoidal dissipation, and the model coefficients are a 2 = 0.4, a 3 = 0.2. 
We are now in the process of applying the pressure-dilatation model to the compressible shear 
layer and the flat plate boundary layer. 

4 Conclusions 

We have obtained a model for the pressure-dilatation after applying scaling arguments ap- 
plied to a formal solution for homogeneous turbulence. This model has been validated 
by direct comparison with DNS results for pressure-dilatation in homogeneous shear flow 
and isotropic turbulence. Eq. (18) is a reasonable approximation for inhomogeneous flows 
without walls. However, future refinements may be necessary to capture different physical 
processes of importance in other flows. For example, a process of importance in shock- 
turbulence interaction and engine flows is a compressive mean velocity field. In such flows, 
the pressure-dilatation model Eq. (18) will have a contribution which is linear in the mean 
compression. Inspection of the dilatation equation for homogeneous compression, obtained 
by taking the divergence of the momentum equation, shows that there should be such a 
term in the pressure-dilatation. Of course, even though the form of the dependence on mean 
compression is already present in our model, in order to get the coefficient of the dependence 
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right, it may be necessary to to add a term like a 4 pKu mim M ] to the model, with a 4 perhaps 
being calibrated with respect to DNS of homogeneous compression. In wall-bounded flows, 
the pressure-dilatation is probably smaller than the estimate given by Eq. (18) because, first, 
as shown by Kim and Lee 10 the rms rapid pressure (except very near the wall) is a smaller 
fraction of ~pq 2 in wall bounded flows compared to homogeneous shear flow, and second, the 
velocity component normal to the wall is preferentially damped. 

The model presented here for the pressure-dilatation is currently being applied in the 
calculation of inhomogeneous flows. As pointed out in the previous paragraph, some exten- 
sions may be necessary in order to incorporate different physical processes in other flows. 
We plan to further develop the pressure-dilatation model for more general flows. 
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Table 1: Parameters for the DNS cases of homogeneous shear flow 


Case 

5 

EM 

Mtfi 

R\, o 

P' 

d' 

mm 

10 


0.3 

24 

0 

0 


15 

gum 


24 

0 

0 

S3 

15 

■rjpi 




0 


Table 2: Parameters for the DNS cases of decaying isotropic turbulence 


Case 

V 


EH 

m 

d! 

D1 

11^11 

E9 

EH 


0 

D2 


1 

mm 

0 

0 

D3 

EM 


27 

0 

0 
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Fig. 2a. Comparison of jrd! and p r d! for Case SI. 
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Fig. 2b. Comparison of the time integrals of p^cL' and p r d' for Case SI. 
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Fig. 3a. Comparison of p c dl and p d! for Case S2. 



Fig. 3b. Comparison of the time integrals of p dl and p d! for Case S2. 
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Fig. 4a. Comparison of p d! and p d 1 for Case S3. 



Fig. 4b. Comparison of the time integrals of p C d' and p r d' for Case S3. 
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Incompressible pressure— dilatation in decaying isotropic turbulence. 



Fig. 7b. Incompressible pressure— dilatation in homogeneous shear turbulence. 
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